From the bulk to monatomic wires: An ab-initio study of magnetism in Co systems 

with various dimensionality 
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A systematic ab-initio study within the framework of the local-spin-density approximation in- 
cluding spin-orbit coupling and an orbital-polarization term is performed for the spin and orbital 
moments and for the X-ray magnetic circular dichroism (XMCD) spectra in hep Co, in a Pt sup- 
ported and a free standing Co monolayer, and in a Pt supported and a free standing monatomic Co 
wire. When including the orbital-polarization term, the orbital moments increase drastically when 
going to lower dimensionality, and there is an increasing asymmetry between the L2 and L3 XMCD 
signal. It is shown that spin and orbital moments can be obtained with good accuracy from the 
XMCD spectra via the sum rules. The (Tz) term of the spin sum rule is surprisingly small for the 
wires, and the reason for this is discussed. 
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The modern methods to prepare nanostructured sys- 
tems made it possible to investigate the influence of di- 
mensionality on the magnetic properties. A central ques- 
tion thereby is how the qualitative behavior will change 
when going from two-dimensional to one-dimensional sys- 
tems because it has been predicted that there is no long- 
range magnetic order at finite temperature in infinitely 
extended one-dimensional systems with short-range mag- 
netic interactions. In the past, several experimental in- 
vestigationSjOf monolayer nanostripes of Fe on vicinal sur- 
faces of WtrB or CuEl have been performed, with a stripe 
width, down to 1-10 nm. Most recently, Gambardella 
et alu succeeded to prepare a high density of parallel 
atomic chains along steps by growing Co on a high purity 
Pt(997) vicinal surface in a narrow temperature range. 
The magnetism of the Co wires was investigatedQ by 
the X-ray magnetic circular dichroism (XMCD). Below 
a blocking temperature a long-range magnetic ordering 
owing to the presence of anisotropy barriers was found 
on the time scale of the experiment. Applying a simplp 
model of exchange coupled superparamagnetic clustersO, 
the anisotropy energy could be obtained from the shape 
of the magnetization curve above the blocking tempera- 
ture and it appeared to be much larger than the one for 
a Co monolayer on Pt which — in turn — is much larger 
than the one of hep Co. Accordingly, a large orbital mo- 
ment of 0.68 /iB per Co atom was found, the highest value 
ever reported for a 3d itinerant electron system. 

So far, magnetism in quasi-one-dimensional systems 
was studied mainly in the insulating material CsNiFa 
where the magnetic Ni ions are arranged along linear 
chains which are well separated from each other so that 
the interchain interaction is only 10"'^ or less of the intra- 
chain interaction. Because of the one-dimensional char- 
acter of this spin system and the easy-plane anisotropy, 
magnetic solitons play an important role for the dynam- 
ical and thermodynamical behavior which h^ been in- 
vestigated by neutron scattering experimentsQ. The dis- 



covery of magnetism in one-dimensional metallic systems 
opens up the chance to extend the research in many 
respects. First, by considering various vicinal surfaces 
the distance between the steps and hence the chains can 
be modified so that the transition from the well-isolated 
chains to interacting chains can be studied. Second, it is 
possible to grow, e.g., biatomic wires along the steps and 
to manipulate the wire length, and to study the respec- 
tive influence on magnetic properties. Finally, we expect 
that the damping of magnetic excitations is larger for 
metals than for insulators and especially large for one- 
dimensional metals, possibly leading—to peculiar proper- 
ties of the nonlinear spin excitationsu. 

The already existing ab-initio calculations for 
monatonMe-|tM|ftsition metal wires considered magnetic 
momentdallOllMBEJ, exchange couplingal3 and magnetic 
anisotropiestllE3liJ. It turned outt3 that in these 
one-dimensional systems the orbital correlations were 
essential and that the magnetic anisotropy and hence 
also the orbital moments could not be calculated reliably 
by the local-spin-density approximation (LSDA) in 
combination with spin-orbit (SO) coupling. Instead, a 
term correcting explicitly for the orbital correlations 
has to be added. It is one of the objectives of the 
present paper to work out the growing importance of the 
orbital correlations with decreasing dimensionality bjj 
calculating orbital (and spin) moments by the LSDAEj 
including SO coupling (in a self-consistent jManner) 
without and with the orbital-polarization termtj which 
takes into account at least in part the orbital-correlation 
effects. The calculations are performed for Co atoms 
in hep Co, in a Pt supported and in a free standing 
monolayer, and in a Pt supported and in a free standing 
monatomic wire. Another objective is to calculate the 
respective XMCD spectra and discuss the influence of 
dimensionality on the accuracy of the spin and orbital 
moments when these ar^jpebtained from the XMCD 
spectra via the sum rulesll3li3. 
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with the XMCD signal = — and with fit = 
fi^ + + fi°. Here A'^h is the number of d holes and (T^) 
is the expectation value of the magnetic dipolar operator 



1 



[a - 3f (f ■ a-)]. 
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FIG. 1: The supercell used to model the monatomic Co wire 
on a vicinal Pt(997) surface. 



To focus on the pure effect of the dimensionality we fix 
the nearest-neighbor distance of the atoms for all consid- 
ered systems to the one of fee Pt (2.77A). For the mono- 
layers and wires we perform supercell calculations, i.e., 
large unit cells which model the structure are repeated 
periodically. In the case of the monolayer the supercell 
consists of two Pt and one Co {111} layers in the fee 
stacking and a vacuum sheet corresponding to two empty 
layers. The vicinal Pt(997) surface with the Co wire at 
the steps is modeled by the supercell shown in Fig. 1 with 
two additional vacuum layers on the top. The supercell of 
the free standing wire is constructed by removing all Pt 
atoms from this supercell. Test calculations have shown 
that the results change only very slightly when going to 
larger supercells. For the monolayers and the wires we as- 
sume perpendicular magnetization, for hep Co it is along 
the c axis. 

We use the tight-binding linear-muffin-tin-orbital 
(LMTO) method in the atomic sphere approximatiorO 
in which we have implemented the SO coupling and 
the tools for calcwJating the XMCD spectra, and 
the WIEN97 codeB which adopts the full-potential 
linearized-augmented-plane-wave (FLAPW) mcthodEll in 
which the SO coupling and the tools for the calculation of 
the magnetooptical effects and XMCD spectraEjS and 
the orbital-polarization termcJ have been implemented. 
The magnetic orbital moments m; — — /iB(^z) and spin 
moments nig = —fiBicz) are calculated direct]^ from the 
wave functions as well as via the sum rulesEjilj from the 
absorption coefficients A^^(e), ^~(e) and n°{e) for hght 
with right-circular, left-circular and z-axis polarization 
at the L2 and L3 edge according to 



where cr denotes the vector of the Pauli matrices. The 
quantities Ep and Ec denote the Fermi energy and a 
cutoff energy. For details for such type of calculations 
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see Ref. |25| The (T^) term is negligible for cubic sur- 
roundings but it is expected to become more and more 
important when reducing the dimensionality of the sys- 
tem. Since it is very difficult to measure (T^), this term is 
often neglected in the spin sum rule when analyzing the 
experimental data. One of the objectives of this paper is 
to asses this critically. 

The results of the LSDA calculations with SO coupling 
are given in the upper part of Table 1. It should be noted 
again that thereby we use the same nearest-neighbor dis- 
tance between the atoms, i.e., the one of fee Pt, for all 
structures. As explained in Ref. the spin and or- 
bital moments derived experimentally via the sum rules 
correspond essentially to the respective part of the va- 
lence wavefunctions with 3d angular momentum, if the 
"background" due to all additional contributions to the 
experimental absorption spectra is subtracted appropri- 
ately. We therefore compare in Table 1 the moments as 
calculated directly from the 3d spin and orbital densities 
with those obtained from the sum rules when taking into 
account for the absorption spectra and for the (Tj) term 
also only that part of the valence wave functions which 
has 3d character. For comparison, we give also the val- 
ues of directly calculated moments including all angular 
momentum contributions. 

For all the structures there is only a slight difference 
between the directly calculated moments with only d con- 
tributions and with all contributions. The spin moments 
increase only moderately with decreasing dimensionality. 
For the wire on Pt our directly calculated FLAPW s pin 
value of 2.06 /ib agrees well with the one of Ref. |l^. 
There is also a satisfactory agreement between the di- 
rectly calculated spin moments and those obtained from 
the spin sum rule when including the (T^) term. The 
{Tz) contribution is negligible for hep Co. In contrast, 
for the free and the Pt supported monolayers, the val- 
ues of ((Tz) obtained from the sum rule when neglecting 
(Tz) are between 20% and 30% smaller than the directly 
calculated values. Astonishingly enough, when going to 
even stronger reduced dimensionality, i.e., for the Pt sup- 
ported wire and even more for the free wire, the (Tz) term 
again appears to be of minor importance. To figure out 
the reason for this startling result. Fig. 2 shows the en- 
ergy distributionc3 for {Tz) for all the structures, i.e., the 
integral over the energy resolved contribution to (Tz) up 
to an energy E (for the integrals up to the Fermi level, 
E — Ep, (Tz) equals the quantity in eq. (2) ). It be- 
comes obvious that for the free standing monolayer and 
for the free standing wire, (Tz) depends drastically on 
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TABLE L The spin and orbital moments (in Bohr magnetons) as calculated from the LSDA without the orbital-polarization 
term (upper part) and with the orbital-polarization term (lower part). For the meaning of the various quantities, see text. The 
upper part gives also the (Tz) term appearing in the spin sum rule. 
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FIG. 2: The energy distribution of (Tz) for hep Co (thin 
dotted line), a Co monolayer on Pt (thin dashed line), a free 
standing Co layer (thick dashed line), a monatomic Co wire 
on Pt (thin solid line) and a free standing monatomic Co wire 
(thick solid line). 



E resulting in sharp structures in the energy distribu- 
tion due to the two- or one-dimensional character of the 
systems. Going to the Pt supported monolayer and wire, 
the energy dependence of (T^) is less pronounced and the 
sharp structures are smeared out, indicating that these 
systems are not really two- or one-dimensional due to the 
presence of the substrate. However, (T^) obtained from 
the integral up to is very small for the free standing 
Co wire but larger for the other systems (except for hep 
Co). This clearly demonstrates that the expected rule 
of thumb — a larger {Tz) for the smaller dimensional- 
ity — does not necessarily hold. We have performed the 
same calculation for a free standing Ni wire. The energy 
distribution of (T^) is very similar but the Fermi level 



is shifted and (Tz) is consequently much larger than for 
the Co wire. This again demonstrates that the symmetry 
arguments alone are not able to estimate the size of {Tz). 

The orbital moments obtained by the LSDA calcula- 
tion with SO coupling are very similar (about O.lS/ie) for 
hep Co and for the Pt supported Co monolayer and wire, 
while the values are sHghtly larger (about 0.18/xb) for the 
free monolayer and considerably larger for the free wire. 
Thereby, there is a good agreement between the directly 
calculated orbital moments and those obtained from the 
calculated spectra via the sum rules. The rather small 
values for the Pt supported monolayer (0.13 /xb) and the 
Pt supported wire (0.14/iB) are iru;onflict with the large 
corresponding experimental valuesB (0.31/iB and 0.68/j,b). 
As discussed in Ref. |2[ the reason is presumably that 
LSDA -I- SO coupling does not appropriately account for 
orbital correlations. We therefore have redone all the cal- 
culations by including in addition the orbital-polarization 
termta. Table 1, lower part, gives the so obtained results 
for the orbital moments (the spin moments are nearly 
unaffected by this additional term). It should be noted 
that the calculations for hep Co have been performed for 
the nearest neighbor distance of fee Pt which explains 
the large moment of 0.4 /ib. For the Co monolayer on 
Pt (Co wire on Pt) our value 0.43 up (0.92 /ib) is even 
larger than the experimental valueO 0.31 /ib (0.68 /xb) 
which may be in part due to the fact that our calcula- 
tions did not take into account any structural relaxations 
at the surface. For the free standing wire, the orbital mo- 
ment appears to be extraordinarily high (2.3 /l(b)- Obvi- 
ously the experimental trend — increase of the orbital 
moment with decreasing dimensionality — is well repro- 
duced when taking into account the orbital-polarization 
term. The increasing orbital moments result in an in- 
creasing asymmetry between the calculated L2 and L3 
XMCD spectra (Fig. 3), in qualitative agreement with 
the respective experimental spectra (Fig. 2 of Ref. ||). 

Both, for the LSDA calculations with and without the 
orbital-polarization term, the directly calculated orbital 
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FIG. 3; The calculated XMCD spectra for hep Co (upper 
part), Co monolayer on Pt (middle part) and Pt supported 
Co wire (lower part). The dashed (solid) lines show the result 
from the LSDA without (with) the orbital-polarization term. 



moments agree very well with those obtained from the 



sum rule. Therefore, for Co in structures of various di- 
mensionality the spin and orbital moments probably may 
be obtained with high accuracy from the sum rules (if in- 
cluding the {Tz) term in the spin sum rule). 

To conclude, we have shown that the orbital corre- 
lations have to be taken into account explicitly in or- 
der to reproduce the experimentally observed trend to 
higher orbital moments and larger asymmetries between 
the L3 and L2 XMCD spectra in Co systems with re- 
duced dimensionality, and we are about to do this also 
for a systematic study of the magnetic anisotropy. In 
addition, we calculate ah initio the further parameters 
appearing in a model Hamiltonian for the monatomic 
wires, i.e., exchange interactions from the adiabatic spin 
wave spectraEj, and the Gilbert damping factorEa. The 
final objective is to come to a comprehensive description 
of the thermodynamic properties (see Ref. ||) and of the 
damped nonlinear excitations (along the lines of Ref. ||) 
of monatomic magnetic linear chains. 
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